arXiv:1501.03483v2 [cond-mat.mtrl-sci] 3 Feb 2015 


Origin of the spin reorientation transitions in (Fei_a:COa;)2B alloy^ 

Kirill D. Belashchenko/ Liqin Ke,^ Markus Dane,^ Lorin X. Benedict,^ Tej Nath Lamichhane,^’^ Valentin 
Taufour,^’^ Anton Jesche,^'^ Sergey L. Bud'ko,Paul C. Canfield,and Vladimir P. Antropov^ 

Department of Physics and Astronomy and Nebraska Center for Materials and Nanoscience, 

University of Nebraska-Lineoln, Lincoln, Nebraska 68588, USA 
Ames Laboratory, U.S. Department of Energy, Ames, Iowa 50011, USA 
Lawrence Livermore National Laboratory, Livermore, California 94550, USA 
Department of Physics and Astronomy, Iowa State University, Ames, Iowa 50011, 

USA 

(Dated: 4 February 2015) 

Low-temperature measurements of the magnetocrystalline anisotropy energy K in (Fei_a;Coa;) 2 B alloys are 
reported, and the origin of this anisotropy is elucidated using a first-principles electronic structure analysis. 
The calculated concentration dependence K{x) with a maximum near x = 0.3 and a minimum near x = 0.8 
is in excellent agreement with experiment. This dependence is traced down to spin-orbital selection rules and 
the filling of electronic bands with increasing electronic concentration. At the optimal Co concentration, K 
depends strongly on the tetragonality and doubles under a modest 3% increase of the c/a ratio, suggesting 
that the magnetocrystalline anisotropy can be further enhanced using epitaxial or chemical strain. 


Magnetocrystalline anisotropy (MCA) of a magnetic 
material is one of its key properties for practical applica¬ 
tions, large easy-axis anisotropy being favorable for per¬ 
manent magnets Intelligent search for new materials 
requires understanding of the underlying mechanisms of 
MCA. This can be particularly fruitful for substitutional 
alloys whose properties can be tuned by varying the con¬ 
centrations of their components. The analysis is often 
relatively simple in insulators, where MCA is dominated 
by single-ion terms which can be deduced from crystal- 
field splittings and spin-orbital (SO) selection rules. In 
contrast, in typical metallic alloys the band width sets 
the largest energy scale, and MCA depends on the de¬ 
tails of the electronic structure. 

The (Fei_a:Co 2 ,) 2 B solid solution^i^— (space group 
lA/mcnr^) is a remarkable case in point. Fe 2 B has a 
fairly strong easy-plane MCA, and C 02 B, at low temper¬ 
atures, a small easy-axis MCA. However, the alloy has a 
substantial easy-axis MCA around x = 0.3,— making it 
a potentially useful rare-earth- free^ permanent magnet. 
At a: ~ 0.5 the MCA again turns easy-plane, peaks at 
X = 0.8, and then turns easy-axis close to a; = 1. These 
three spin reorientation transitions must be related to 
the continuous evolution of the electronic structure with 
concentration. The goal of this Letter is to elucidate the 
origin of this rare phenomenon. 

First, we report the results of experimental mea¬ 
surements at low temperatures. Single crystals of 
(Fei_a;Coa;) 2 B were grown from solution growth out of 
an excess of (Fe,Co) which was decanted in a centrifugei^ 
The single crystals were grown as tetragonal rods which 
were cut using a wire saw to give them the shape of a rect¬ 
angular prism. The demagnetization factor was calcu- 


“^iThis article has been accepted by Applied Physics 
Letters. After it is published, it will be found at 
http://scitation.alp.org/content/aip/j ournal/apl 


lated using Ref. M Field-dependent magnetization mea¬ 
surements were performed in a Quantum Design MPMS 
at 2 K in fields up to 5.5 T. The MCA energy K was 
determined as the area between the two magnetization 
curves, with the field parallel and perpendicular to the c 
axis, taken at the same temperaturei^ The results shown 
in Fig. [T] (blue stars) measured at 2 K confirm the non¬ 
monotonic concentration dependence, in good agreement 
with the measurements at 77 K from Ref. |3. 



FIG. 1. Calculated (black circles) and experimental (blue 
stars) MCA energy K in (Fei-a;Coa;) 2 B alloys. Gray curve: 
Kso- The other lines show the spin decomposition K„„t. 


Density-functional calculations using several different 
methods show that the choice of the exchange-correlation 
potential and other computational details strongly affect 
the calculated MCA in Fe 2 B and C 02 B. We have ascer¬ 
tained that this sensitivity is largely due to the variation 
of the exchange splitting, which controls the position of 
the Fermi level relative to the minority-spin bands. The 
systematic variation of MCA with x is also controlled 
by the Fermi level shift. This continuous variation is, 
therefore, reliably predicted as long as the end points are 
correctly fixed using experimental input. 
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The results reported below were obtained using the 
Green’s function-based formulation of the tight-binding 
linear muffin-tin orbital (GF-LMTO) method.— Substi¬ 
tutional disorder was treated using our implementation 
of the coherent potential approximation (GPA),— with 
the SO coupling included perturbatively as described in 
the supplementary material.— 

The lattice constants and the internal coordinates are 
linearly interpolated between the experimental data for 
the end compounds extrapolated to zero temperatures^ 
a = 5.109 and 4.997 A, c = 4.249 and 4.213 A, and 
u = 0.166 and 0.168 for Fe 2 B and C 02 B, respectively. 
The exchange and correlation are treated within the gen¬ 
eralized gradient approximation (GGA)ji^ 

The correct exchange splitting at the end points can 
be enforced by using the experimental data^i — for the 
magnetization M: 1.9 /rs/Fe in Fe 2 B and 0.76 fis/Co 
in G 02 B. In Fe 2 B it is only slightly overestimated, but 
in C 02 B it is much too large at about 1.1 hb/Co in all 
density-functional calculations. The relatively small spin 
moment of Co indicates a pronounced itinerant charac¬ 
ter of magnetism in C 02 B, which tends to be sensitive to 
quantum spin fluctuations. Therefore, we introduced a 
scaling factor for the local part of the effective magnetic 
field in the GGA functional for the Co atoms and ad¬ 
justed it to match the experimental magnetization. The 
resulting scaling factor of 0.80 was then used for Co atoms 
at all concentrations. The spin moments on different 
atoms and the total spin magnetization obtained in this 
way are shown in Fig. [2l While the Fe spin moment is 
almost constant, the Co spin moment declines with x. 
Moreover, this decline accelerates at x ^ 0.6, which is in 
excellent agreement with experimental data.— 
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FIG. 2. Spin moments on different atoms and the total spin 
magnetization M per transition-metal atom. 

The MCA energy K was obtained by calculating the 
single-particle energy difference for in-plane and out-of- 
plane orientations of the magnetization while keeping the 
LMTO charges fixed at their self-consistent values found 
without SO coupling. A uniform mesh of 30^ points in 
the full Brillouin zone provided sufficient accuracy for the 
k integration. We have verified that the values of K for 


pure Fe 2 B and C 02 B agree very well with Hamiltonian 
LMTO results. The concentration dependence of K is 
shown in Fig. [TJ The agreement with low-temperature 
experimental data is remarkably good, suggesting that 
the electronic mechanisms of MCA are correctly cap¬ 
tured in the calculations. If the spin moment of Co is 
not corrected by scaling the exchange-correlation field, 
the downward trend in K at the Co-rich end continues 
to large negative values in disagreement with experiment. 

In 3d systems the SO band shifts are usually well 
described by second-order perturbation theory, except 
perhaps in small regions of the Brillouin zone. Conse¬ 
quently, when MCA appears in second order in SO cou¬ 
pling (as in the tetragonal system under consideration), 
the anisotropy of the expectation value of the SO opera¬ 
tor AEso = {yso)x — {Vso)z is approximately equal to 
2 AT— (Here x or z shows the orientation of the magne¬ 
tization axis.) We therefore denote Kso = AEso/‘^ and 
use the expression 2(SL) = (Ti')tt - {Lz')u + (-^+)4t + 
to separate the contributions to K by pairs of spin 
channels. Here we use Lz' to denote the component of 
L parallel to the magnetization axis, to avoid confusion 
with the crystallographic 2 ; axis; L± are the usual linear 
combinations of the other two (primed) components of L. 
The contributions to Kso are accumulated as en¬ 
ergy integrals taking into account the energy dependence 
of the SO coupling parameters. The results for Kso and 
ATctct' are shown in Fig. [T] First, we see that Kso is close 
to K, confirming the validity of this analysis. Second, the 
nonmonotonic concentration dependence of K is almost 
entirely due to K^^ for x < 0.7. While K^^ is sizeable, 
it depends weakly on x in this region. Additional discus¬ 
sion about the spin decomposition of MCA is provided 
in the supplementary material— 

Let us now analyze the electronic structure and the de¬ 
tails of SO coupling-induced mixing for the minority-spin 
electrons. To resolve the minority-spin contribution to K 
by wave vector k, we calculated the minority-spin spec¬ 
tral function in the presence of SO coupling and found its 
first energy moment at each k. Fig. [3] shows the differ¬ 
ence of these integrals for magnetization along the x and 
2 ; axes at three key concentrations: pure Fe 2 B {x = 0), 
the maximum of K {x = 0.3) and its minimum (x = 0.8). 
We have checked that the Brillouin zone integral of the 
k-resolved MCA energy (summed up over both spins) 
agrees almost exactly with the value of K calculated in 
the usual way. 

Fig. [3] shows that the MCA energy accumulates over a 
fairly large part of the Brillouin zone. At x = 0 negative 
contributions to K dominate over most of the Brillouin 
zone. At X = 0.3 both positive and negative contribu¬ 
tions are small. At x = 0.8 there are regions with large 
positive and large negative contribntions. Overall, it ap¬ 
pears that the most important contributions come from 
the vicinity of the FXM {kz = 0) plane and from the 
vicinity of the FH [k^ = ky = 0) line. 

The partial minority-spin spectral functions for the 
transition-metal site are displayed in Fig. 0] (panels (a)- 
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FIG. 3. Brillouin zone map of the k-resolved minority-spin contribution to K at: (a) x = 0, (b) x = 0.3, (c) x = 0.8. Half of 
the Brillonin zone is shown; the top face of the plot is = 0. Points H (same as M) and X are shown in panel (b). The color 
intensity indicates the magnitnde of negative (blue) and positive (red) valnes. 


(c)) along the important high-symmetry directions for 
the same three concentrations. The coloring in this fig¬ 
ure resolves the contributions from different 3d orbitals. 
At a; = 0 the spectral function resolves the conventional 
electronic bands of pure Fe 2 B (an imaginary part of 0.004 
eV is added to the energy to acquire them). At x = 0.3 
and X = 0.8 substitutional disorder broadens the bands 
by a few tenths of an electronvolt, but their identity is 
in most cases preserved. Thus, we will discuss the SO- 
induced band mixing in the alloy, bearing in mind that 
band broadening should reduce the values of MCA at 
intermediate x. 

As we have learned above, the dominant concentration 
dependence of K comes from the term, where z' 

is the magnetization axis. The electronic states on the 
whole rXM plane can be classified as even or odd under 
reflection z —>■ —z. Even (odd) states have m = 0, ±2 
(m = ± 1 ) character and appear red and green (blue) in 
Fig. m States of different parity do not intermix on this 
plane in the absence of SO coupling, as is clearly seen 
in Fig. 0^. The selection rules for SO coupling of the 
minority-spin states follow from the definite parity of the 
components of L under reflection. (even; relevant for 
M II z) only mixes states of the same parity on the FXM 
plane, or more generally orbitals of the same m. Coupling 
between states of the m = ±2 character (red) is stronger 
compared to states of the m = ±1 character (blue). In 
contrast, (odd; relevant for M || x) couples states of 
the opposite parity on the FXM plane, or more generally 
orbitals m and mil. All these couplings contribute to 
K only when the Fermi level lies between the two states 
that are being coupled. Whenever or couples such 
states, there is a negative contribution to the energy of 
the system with the corresponding direction of M. 

With the help of Fig. 0]we can now deduce which cou¬ 
plings contribute to K at different concentrations. At 
X = 0 the Fermi level lies between the filled even states 
(bands 1-2) and empty odd states (bands 3-4) near F. 
Coupling of these states by contributes to negative 
K. Filling of the hole pocket at F (bands 3-4) with in¬ 
creasing X suppresses this contribution. At a: = 0.3 the 
odd bands 3-4 are filled (Fig.0]D), and the minority-spin 


contribution to K is small (Fig.[5]D). These two cases are 
sketched in Fig. |4ji. 

At still larger x the odd band 5 gets gradually filled, ac¬ 
tivating the negative contribution to K from the mixing 
of band 5 with empty even bands 6-7. This trend contin¬ 
ues till about X = 0.8, where an even pair of bands 6-7 
(degenerate at F) begins to fill (Fig. 0):). Mixing of these 
bands by Lz leads to an intense positive contribution to 
K near the F point (clearly seen in Fig. [3]:), and the trend 
reverses again. Thus, the nonmonotonic dependence of 
K in the whole concentration range is explained. Note 
that if the exchange-correlation field on the Co atoms 
is not scaled down to bring the magnetization in agree¬ 
ment with experiment, the exchange splitting remains 
too large, and bands 6-7 remain empty up to cc = 1. As a 
result, without this correction the uninterrupted negative 
trend brings K to large negative values in disagreement 
with experiment. 

To assess the effect of atomic relaxations on K, we op¬ 
timized all inequivalent structures with two formula units 
per unit cell using the VASP codeia and the GCA. The 
volumes were fixed at the same values as in the CPA 
studies at the same x, while the cell shape and inter¬ 
nal coordinates were relaxed. Since all these supercells 
preserve the az reflection plane, the displacements of Fe 
and Co atoms are confined to the xy plane. All displace¬ 
ments were less than 0.025 A. One of the two structures 
at X = 0.5 breaks the Ci symmetry. For this structure 
we took the average of the energies for M || x and M || y 
as the in-plane value in the calculation of K, which cor¬ 
responds to the averaging over different orientations of 
the same local ordering. 

The changes in the absolute value of K due to the re¬ 
laxation and its values (meV/f.u.) in the relaxed struc¬ 
tures were found to be: —26% and —O.II in Fe 2 B, — 6 % 
and 0.25 in Fe 1 . 5 Coo. 5 B, —13% and 0.15 in the FeCoB 
[100] superlattice, 6 % and 0.12 in the FeCoB [110] su¬ 
perlattice, — 11% and —0.31 in Feo. 5 C 01 . 5 B), and —19% 
and —0.04 in C 02 B. While MCA tends to be larger for 
ordered structures, the concentration trend in supercell 
calculations agrees well with the CPA results for disor¬ 
dered alloys. Although this set of unit cells is limited, the 
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FIG. 4. (a-c) Minority-spin partial spectral functions for the transition-metal site in the absence of SO coupling at (a) a; = 0, 
(b) X = 0.3, and (c) x = 0.8. Energy is in eV. (d) Level diagram and SO selection rules at the F point (bands 1-4). Color 
encodes the orbital character of the states. The intensities of the red, blue and green color channels are proportional to the 
sum of xy and x^ — {m = ±2), sum of xz and yz (m = ±1), and z^ (m = 0) character, respectively. 


results suggest that local relaxations do not qualitatively 
change the concentration dependence shown in Fig. [T] 
Larger positive values of K are favorable for perma¬ 
nent magnet applications. Our electronic structure anal¬ 
ysis shows that the maximum near x ~ 0.3 corresponds 
to the optimal band filling. Further raising of K re¬ 
quires favorable changes in the band structure, which 
could be induced by epitaxial or chemical strain. We 
therefore considered the dependence of K a,t x = 0.3 on 
the volume-conserving tetragonal distortion. The results 
plotted in Fig. [S]show a very strong effect: K is doubled 
under a modest 3% increase in c/a due to the sharply 
increasing spin-flip contributions. A more detailed anal¬ 
ysis shows that the latter is largely due to the increase in 
the c parameter. On the other hand, the minority-spin 
contribution increases with decreasing volume. Thus, in¬ 
creasing c and decreasing a both have a positive effect on 
MCA. This enhancement could be achieved through epi¬ 
taxial multilayer engineering. The search for a suitable 
alloying element to enhance the c/a ratio is an interesting 
subject for further investigation. 



FIG. 5. MGA energy A as a function of the c/a ratio at 
X = 0.3 (the value c/a = 1 is assigned to the unstrained 
lattice). Kso and Ao-o-' are also shown. 

In conclusion, we have explained how the spin re¬ 
orientation transitions in (Fei_a;Coa;) 2 B alloys originate 
in the SO selection rules and the consecutive filling of 


minority-spin electronic bands of particular orbital char¬ 
acter. Near the optimal 30% concentration of Co, the 
MCA energy is predicted to increase quickly with the c/a 
ratio, which could be implemented by epitaxial strain or 
a suitable chemical doping. 
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Origin of the spin reorientation transitions in (Fei_^Coa;)2B alloys 
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I. IMPLEMENTATION OF SPIN-ORBIT COUPLING IN GREEN’S FUNCTION-BASED LINEAR MUFFIN-TIN 
ORBITAL (LMTO) METHOD 

Magnetocrystalline anisotropy (MCA) of substitutional alloys is often studied using the fully-relativistic Korringa- 
Kohn-Rostocker or Green’s function-based LMTO (GF-LMTO) method combined with the coherent potential approx¬ 
imation (CPA)fii“— Here we employed a perturbative implementation of spin-orbit coupling in GF-LMTO. Dealing 
with conventional spinor wave functions, it simplifies the analysis of MCA without a significant loss of accuracy 
compared to the fully-relativistic approach. 


A. Relation between the Hamiltonian and Green’s function-based LMTO formulation 


The so-called orthogonal Hamiltonian and Green’s function-based LMTO formulations are equivalent to second 
order in {E — E,j) where E^ is the linearization energyThe Hamiltonian in the orthogonal basis is 

Horth = C + ^/Xs{l--tS)-^'/K, ( 1 ) 

where C, A, and 7 are the diagonal matrices containing the LMTO potential parameters, and S is the structure 
constant matrix. On the other hand, in GF-LMTO the Green’s function 

G{z) = \{z) E ^x{z)[P{z) - S']"VW, (2) 


is constructed using the diagonal matrices containing the potential functions 

z-C , , TA 


P{z) = 


fl{z) = 


A(z) = 


+ 7i(z — G)’ A + 7(z — C) ’ A -I- 7(z — G) 

These definitions are related through energy derivatives (denoted by an overdot): 

M^ = P, A = -ip/A 


(3) 

(4) 


These relations guarantee that G(z) does not have poles at the points where P(z) has simple poles. It is straightforward 
to show that G(z) = (z — fforth)~^, which means that G(z) in ([2]) is the exact resolvent of Horth- 

Third-order parametrization of the potential functions^ is often considerably more accurate compared to the second- 
order one. The accuracy of this parametrization is similar to that of the three-center LMTO Hamiltonian, although 
they are not exactly equivalent. In the nearly orthogonal LMTO basis, the Hamiltonian and overlap matrices are^ 


P^ — P^orth ~t“ (Flori/i Pu^PuP(^P^orth A^/), 

0 = 1 + (Plorth — Ev)p{Plorth ~ Ev) (5) 

where Ei, and p are diagonal matrices containing the linearization energies and the LMTO parameters p = {4^), and 
Horth is the same as in ©• The generalized eigenvalue problem is then written as 

det [{Horth - E)- {Horth - E^){E - E„)p{Horth - E^)] = 0 (6) 

The third-order term can be treated as a perturbation, and the variational correction formally amounts to the 
substitution El ^ El + {E — E^)^p in the second-order eigenvalue equation det{Horth — E) = 0. 

In the GF-LMTO formulation, third-order accuracy is similarly achieved by redefining the potential functions)^!^ 

P = P{z), il = y/pp.{2), A =-bz'A(z) (7) 

where z = z+p{z — E^)'^ and z' = dz/dz. These definitions are designed to preserve the relations (|3]). If the third-order 
Green’s function G is defined similarly to (I2|) but with the third-order potential functions ([7]), it can be shown that 

G = -\^ + VP{z - Horth)~^'/I'. 
z z 


( 8 ) 
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A solution of det{Horth — z) = 0 corresponds to a pole of G, and the factors \/|' guarantee that the residue at that 
pole is equal to 1. Thus, the eigenstates of the three-center LMTO Hamiltonian are correctly represented by G to 
third order in E — E^. The first term in ([8|) introduces two unphysical poles at z = E^, ±ify/^ for each orbital. It is 
a diagonal matrix which is real on the real axis and has no effect on the spectrum of the physical states. 


B. Spin-orbit coupling in GF-LMTO 

In the Green’s function formalism the perturbation can be introduced through the Dyson equation 

Gv=Go + Go^Gv (9) 

where Gq and Gv are the unperturbed and perturbed Green’s functions and E the energy-dependent self-energy due 
to the perturbation V. For V representing the spin-orbit coupling, E is local. 

In the second-order representation Go = {z — Horth)~^, and therefore Gy = {z — Horth — E)“^. This Gy can be 
constructed from the potential functions by denoting 

Cy = C + i: (10) 


and defining the following nondiagonal matrices: 


Py = y/E[A + {z- Gy)-f]-\z - Gy) 


Va 


"fE'/A{z-Cy) 


•n -1 


/i^ = [A + 7 (z - Cy)]-iyA = {z- Gy)-^srAPy 
= yA[A + (z - Gy)-l]-^ = PyVA{z - Gy)-^ 
Ay = [A -F 7(z - C'y)]"^7 = 7[A + {z- Gy)-i]~^ 


The relations (0]) are then replaced by 

Pv = MyMyj Ay = --{py) ^P{Pv) 
A tedious but straightforward derivation then shows 

Gy = Ay + My(^y-5)-Vy- 


( 11 ) 


( 12 ) 


(13) 


If the unperturbed Green’s function is calculated to third-order accuracy, we can still proceed from Dyson’s equation, 
but the situation is complicated by the first term k = —z"j{2P) in Eq. (|8]) for Gq, which modifies the perturbed 
Green’s function Gy. However, for spin-orbit coupling this modification may be neglected, as we now argue. Let us 
denote Go = k -I- Gq and write down the diagrammatic expansion of Gy in powers of E. It is just a usual expansion 
in which each Green’s function line can be either Go or k. Now let us resum all diagram insertions with no external 
lines and no Go lines inside. This gives a renormalized self-energy E = (1 — Ek)“^E. Taking this into account and 
performing the remaining summations, we find 

Gy = G -I- K -I- kEk -I- kSG -I- GSk (14) 


where 


G= (l-GoE)-iGo. 


(15) 


The quantity Ek ~ p{E — Ei,)T, is very small for spin-orbit coupling in transition metals: E ~ 50 meV, p ^ 0.01 
eV“^, E — E,j ^ 5 eV. Since k is analytic on the real axis, E has no poles there too. Therefore, the poles of Gy 
coincide with the poles of G. The latter are just the poles of Go which are shifted by the self-energy E. But E differs 
from E by a factor (1 — Ek)“^ which is equal to 1 up to small corrections of order Ek. Up to similar corrections, the 
residues of the poles of Gy are equal to those for G. Thus, neglecting all small terms or order Ek compared to 1, we 
find Gy ~ G -I- K. 

Substituting the second term in ([8]) as Go in dig, we find the perturbed Green’s function: 


Gy = 



Horth 



'/P. 


(16) 
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where we have dropped the inconsequential term k. This expression for Gy is equivalent to 

Gv = Av + Jly{Pv — S) ^fly (17) 

with redefined potential functions 

py = Py{z), jly ='/i'y,y{z), fly = ^y{z)y/I', \y = \y {z)\fp, (18) 

where instead of Gy in m we should use 

Cy = G + (19) 

The self-energy matrix for spin-orbit coupling within the given I subspace is defined as 

'P'SO = ii<Tcr'{z){lma\Sli\lm'a') (20) 

where the energy dependence of the coupling parameter ^lacr' {z) comes from the radial wave functions. To second 
order in (z — the radial function in the LMTO basis is 4>uia(,z) = + {z — E^)ip^ia-. Therefore, we have 

^laa'{z) = {(j)ula{z)\Vso{r)\(j)ulc,' {z)) 

- I (^) I ) E {z I (^) I ) E {,Z 

+ iz- E^i„){z- E^i^>){^^i„\Vsoir)\(j)i,l^>)]/y/[l +Pla{z - EyiaY][l + pia’{z - Eyi^'Y] (21) 

where the square roots in the denominator come from the normalization of ffuiaiz). 

An approximate form can be obtained by neglecting the terms of order p{E — E^p when they appear in the 
perturbation, i. e. approximating •pW r; 1 in (jl9p and neglecting the square root in the denominator of (1211) . We 
found that for (Fei_a;Coa;) 2 B alloys this approximation captures all the essential physics (see, for example. Fig. 1), 
while the p{E — E^Y terms only change K by an average of about 10%. 

The above treatment of spin-orbit coupling is extended to CPA with no modifications. Some implementation notes 
about our CPA code can be found in Ref. d. 


II. RELATION TO SPIN-ORBIT COUPLING ENERGY 

When MCA energy K appears in second order in spin-orbit coupling, the anisotropy of the expectation value of 
the SO operator AEso = {Vso)x — {Vso)z is approximately equal to 2K^ Here we transcribe this relation in terms 
of Green’s functions. 

The single-particle energy is by definition 

Esp = -ImTr J zGydz (22) 

The second-order term from perturbation theory is 

EfJ = -ilmTr J zGo^GoE^Godz = -^ImTr J zG^SGoSdz (23) 

The expectation value of V 501 with the radial integral included in E, in the same order is 

(Vso) = - i ImTr [ T.Gdz « Tr [ SGoEGodz. (24) 

TT J TT J 

An exact quasiparticle Green’s function with simple poles of unit residue satisfies dG/dz = —G^. For Gq this identity 
is satisfied approximately up to a real term of second order in {E — E^- Thus, setting dGo/dz « —Gq, we find 

Gg) [ z~iGo^Go^)dz + {{dE/dz)) = + {{dE/dz)) (25) 

where the first term was integrated by parts, and {{dE/dz)) denotes the terms with the energy derivative of E. For 
substitutional alloys treated in CPA, the Green’s function does not have a purely quasiparticle structure, and in 
general the relation dG/dz = —G^ is violated due to the energy dependence of the coherent potential. This energy 
dependence contributes additional terms to the right side of (1^51) . However, here we are dealing with relatively weak 
substitutional Fe/Co disorder, and the electronic structure largely preserves its quasiparticle character. Therefore, we 
may expect the relation Elf Ri {Vso)/‘^ to hold approximately, which is borne out by calculations. Deviations from 
this relation can also occur near the points of degeneracy close to the Fermi level, where second-order perturbation 
theory breaks down. As seen in Fig. 4, this situation appears near x = 0.8 where the Fermi level cuts through a pair 
of flat, nearly degenerate bands. From Fig. 1 we see the largest difference between K and Kso at this point. 
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III. SPIN-RESOLVED ORBITAL MOMENT ANISOTROPY 


Orbital moments and their dependence on the direction of magnetization are often discussed in connection with 
MCAiS^— In (Fei_a;COa;) 2 B both spin channels contribute comparably to K, and therefore there is no direct relation 
between K and the anisotropy of the total orbital moment. However, the spin-flip terms are small and change 
slowly when x is not close to 1 (see Fig. 1). Since the energy dependence of the SO parameters is relatively weak, 
the contributions and are closely related to the orbital moment anisotropy (OMA) for the states of the 
corresponding spin. Fig. ISII shows that the concentration-weighted average of the spin-down OMA behaves similarly 
to Kii. The spin-up OMA is almost constant except for the Co-rich end where, as seen in Fig. 1, the spin-flip terms 
also increase substantially. Note that negative OMA for spin-up states corresponds to positive but the sign of 
OMA for spin-down states is the same as that of K^. This is because (^ 2 ) 4 , 4 , appears in Kso with a minus sign, as 
reflected in Hund’s third rule. 



FIG. SI. Orbital moment anisotropy (OMA) for majority and minority-spin states of Fe and Co. The thick black line shows 
the concentration-weighted average of Fe and Co contributions from the minority-spin states. 


The concentration dependence of and of the majority-spin OMA is weak because the majority-spin 3d bands 
are fully filled. This feature is typical for Fe, Co, and Ni-based magnetic alloys. Therefore, while K is almost never 
proportional to the total OMA, their derivatives with respect to the concentration may be strongly correlated in many 
magnetic 3d alloys. 
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